###########################CARREGAR PACOTES##############################
#########################################################################
require(RcmdrMisc) #pacote para carregar arquivos do Excel
require(agricolae)
require(dae)
require(ExpDes)


#Mudar diretório
setwd("D:/OneDrive/UFMT/Disciplinas/Veterinaria/experimental/2019/dados")
setwd("C:/Users/ander/OneDrive/UFMT/Disciplinas/Veterinaria/experimental/2019/dados")
source("funcoes.txt")



#########################################################################
#################Análise do DBC em Esquema Fatorial######################
#########################################################################



##Importar dados
dados=readXL("peixe1.xls",sheet = 1,stringsAsFactors=TRUE)
attach(dados)
str(dados)

T=as.factor(temperatura)
F=as.factor(foto)
######Análise de Variância####################################
model=aov(gsi~bloco+T*F)
anova(model)
cv.model(model)

#####PRESUPOSIÇÕES DA ANOVA
#Gerar os erros padronizados do modelo
e=rstandard(model)

#Teste de Kolmogorov-Smirnov
m=mean(e)
dp=sd(e)
ks.test(e,"pnorm",sd=dp)

############Verificar suposição de Homogeneidade de Variância
#teste de Oneill Mathews
Trat=interaction(T,F)
oneillmathewsDBC(gsi,Trat,bloco)

##Teste de Durbin-Watson - Independencia
set.seed(12345)
durbinWatsonTest(model,reps=10000)

fat2.crd(Levorfanol,Epinefrina,VR1,fac.names = c("Levorfanol","Epinefrina"),quali=c(TRUE,TRUE),mcomp = "tukey")



###Teste de aditividade dos efeitos de tratamento e bloco
df<-df.residual(model)
MSerror<-deviance(model)/df
nonadditivity(gsi, bloco, Trat, df, MSerror)


fat2.rbd(foto,temperatura,bloco,gsi,fac.names = c("foto","temperatura"), quali = c(TRUE,FALSE),mcomp = "sk")
tapply(gsi,list(T,foto),mean)

grafico.regfat2(foto,temperatura,gsi,grau=c(3,3),ylab="GSI",xlab="Temperatura",nomes=NULL,local="topright")